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I investigate numerically the phase transitions of two-component generalizations of binary spread- 
ing processes in one dimension. In these models pair annihilation: AA — > 0, BB —* 0, explicit particle 
diffusion and binary pair production processes compete with each other. Several versions with spa- 
tially different productions have been explored and shown that for the cases: 2A — > 3A, 2B — > 3B 
and 2A — > 2AB, 2B — > 2BA a phase transition occurs at zero production rate (a = 0), that belongs 
to the class of N-component, asymmetric branching and annihilating random walks, characterized 
by the order parameter exponent (5 = 2. In the model with particle production: AB — > ABA, 
BA — * BAB a phase transition point can be located at <r c = 0.3253 that belongs to the class of the 
one-component binary spreading processes. 



One-dimensional, non-equilibrium phase transitions 
have been found to belong to a few universality classes, 
the most robust of them is the directed percolation (DP) 
class . According to the hypothesis of ||[| all contin- 
uous phase transitions to single absorbing states in ho- 
mogeneous, single component systems with short ranged 
interactions belong to this class provided there is no addi- 
tional symmetry and quenched randomness present. The 
most well known exception from the robust DP class is 
the parity conserving (PC) class g, where a mod 2 con- 
servation of particles happens (example in even offspring 
branching and annihilating random walks (BARWe)) and 
in multi-absorbing state systems where an exact Zi sym- 
metry is satisfied too ||. There are other classes being 
explored recently where the total number of particles is 
conserved |7|-|l^]. 

In multi-component systems bosonic field theory JIJ] 
simulations |l3l and density matrix renormalization 
group analysis|L6| have revealed the universality class of 
the generalization of the BARWe class. Hard-core parti- 
cle exclusion effects can change both the dynamic [ |n]Jl8| 
and static p^ , pj| -p2]| behavior of one dimensional sys- 
tems by introducing blockades into the particle dynamics. 
Earlier it was shown that infinite number of conservation 
laws emerge in stochastic deposition-evaporation models 
of Q-mers in one dimension (2^JU] that split up the phase 
space into kinetically disconnected sectors. That results 
in initial condition dependent autocorrelation functions. 

In [jl8| the two-component generalization of the anni- 
hilating random walk model was introduced taking into 
account hard-core repulsion of particles (2-ARW): 

AA^% BB^% A%A%A B%A$B AB ^ BA 

(1) 

(where A and d denote the annihilation and diffusion 
rates) and showed that the initial conditions influence 
the particle density (order parameter) decay and the dy- 
namical exponents. By adding pair creation processes 
(A A- 2BA, B A 2AB) to this model a continuous phase 



transition occurs at a = creation rate and two univer- 
sality classes appear depending on the arrangement of 
the offsprings relative to the parent Namely if the 
parent separates the offsprings: A A BAB (2-BARW2s) 
the steady state density will be higher than in the case 
when they are created on the same site: A A ABB 
(2-BARW2a) for a given branching rate a because in 
the former case they are unable to annihilate with each 
other. This results in different off-critical order param- 
eter exponents for the symmetric and the asymmetric 
cases (P s — 1/2 and j3 a = 2). This is in contrast to the 
widespread beliefs that bosonic field theory can well de- 
scribe reaction-diffusion systems in general. In the field 
theoretical version [[uj , where the AB <-> BA exchange is 
allowed the critical behavior is different. Mean- field like 
and simulation results led Kwon et al. [|l9| to the assump- 
tion that in one-dimensional reaction-diffusion systems a 
series of new universality classes should appear if particle 
exclusion is present. 

In a recent paper |2lJ I showed that if one adds single 
particle creation to the 2-ARW model: 

A-^AB B^BA (2) 

a continuous phase transition occurs again at a = 
and the critical exponents coincide with that of the 2- 
BARW2s model, although the the parity of particle num- 
ber is not conserved. Therefore this conservation law that 
was relevant in case of one component systems (PC versus 
DP class) is irrelevant here. In |2l[ I made a hypothesis 
that in coupled branching and annihilating random walk 
systems of N-types of excluding particles with continuous 
transitions at a = two universality classes exist, those 
of 2-BARW2s and 2-BARW2a, depending on whether the 
reactants can immediately annihilate (i.e. when similar 
particles are not separated by other type(s) of particle (s)) 
or not. These classes differ only by the off-critical expo- 
nents, while the on-critical ones are the same. This is due 
to the fact that the critical point is at zero branching rate 



1 



(a = 0) and therefore they are the ones determined for 
the 2-ARW model @,||. 

In this paper I extend the investigations to coupled 
binary production spreading processes, where new uni- 
versal behavior has recently been reported. Studies on 
the annihilation fission process 2 A — > 0, 2 A — > 3 A, 
A% «-> %A p5|-|29t found evidence that there is a phase 
transition in this model that does not belong to any pre- 
viously known universality classes. This model without 
the single particle diffusion term - the so called pair con- 
tact process (PCP), where pairs of particles can annihi- 
late or create new pairs - was introduced originally by 
Jensen |}2| and while the static exponents were found to 
belong to DP class the spreading ones show non-universal 
behavior. By adding explicit single particle diffusion J2(| 
Carlon et al. introduced the so called PCPD particle 
model. The renormalization group analysis of the corre- 
spondin g b osonic field theory was given by Howard and 
Tauber |2q| . This study predicted a non-DP class tran- 
sition, but it could not tell to which universality class 
this transition really belongs. An explanation based on 
symmetry arguments are still missing but numerical sim- 
ulations suggest [^8],[33| that the behavior of this system 
can be well described (at least for strong diffusion) by 
coupled sub-systems: single particles performing annihi- 
lating random walk coupled to pairs (B) following DP 
process: B — + 2B, B — > 0. The model has two non- 
symmetric absorbing states: one is completely empty, in 
the other a single particle walks randomly. Owing to this 
fluctuating absorbing state this model does not oppose 
the conditions of the DP hypothesis. 

In the low diffusion region (d <~ 0.4) some exponents 
of the PCPD model are close to those of the PC class but 
the order parameter exponent (/3) has been found to be 
very far away from both of the DP and PC class values 
psfl . In fact this system does not exhibit neither a Z2 
symmetry nor a parity conservation that appear in mod- 
els with PC class transition. In the high diffusion region 
the critical exponents seem to be differe nt |2^ , ^8| , |30|] sug- 
gesting an other universality class there [M. This is also 
supported by the pair mean-field results |26|] . A recent 
universal finite size scaling amplitude study ]31| suggests 
however that a single universality class with strong cor- 
rections to scaling may also be possible. 

It is conjectured by Henkel and Hinrichsen |m| that 
this kind of phase transition appears in models where 
(i) solitary particles diffuse, (ii) particle creation requires 
two particles and (iii) particle removal requires at least 
two particles to meet. Very recently Park et al. |55| have 
investigated the parity conserving version of the PCPD 
model (2A -» AA, 2A -» 0, A® «-> %A) and contrary the 
apparent conservation law they have found similar scal- 
ing behavior that led them to the assumption that the 
binary nature of the offspring production is a necessary 
condition for this class. Other conditions that would in- 
fluence the occurrence of this class should be clarified 
too. In this paper I address the question whether the 
particle exclusion effects arc relevant like in the case of 



BARW processes and whether the hypothesis set up for 
N-BARW systems |2l| could be extended. 

One site update step of the applied algorithms consist 
of the following processes. A particle is selected ran- 
domly. A left or right nearest neighbor is chosen with 
probability 0.5. With probability a a pair production is 
attempted in case of an appropriate neighbor. Otherwise 
(with probability d = \ = 1 — 0) a hopping is attempted 
if the neighbor is empty or if the neighbor if filled with 
a particle of same type they are annihilated. The fol- 
lowing models with the same diffusion and annihilation 
terms as ([!]) and different production processes will be 
investigated here. 

a) Production and annihilation random walk model (2- 
PARW): 

AA AAB, AA ^ BAA, (3) 
BB ^ BBA, BB ^ ABB . (4) 

b) Symmetric production and annihilation random 
walk model (2-PARWS): 

AA -2-> AAA, (5) 
BB BBB . (6) 

c) Asymmetric production and annihilation random 
walk model (2-PARWA): 

AB ^ ABB, AB ^ AAB, (7) 
BA BAA, BA ^ BBA . (8) 

d) Asymmetric production and annihilation ran- 
dom walk model with spatially symmetric creation (2- 
PARWAS): 

AB ^ ABA, AB ^ BAB, (9) 
BA ^ BAB, BA ^ ABA . (10) 

The evolution of particle densities were followed by 
Monte-Carlo simulations started from randomly dis- 
tributed A, B, sites in systems of sizes L = 10 5 and 
periodic boundary conditions. 

The 2-PARWA (c) model does not have an active 
steady state. The AA and BB pairs annihilate them- 
selves on contact, while if an A and B particle meets an 
AB — > ABB — > A process reduces out blockades there- 
fore the densities decay with the p oc t~ x / 2 law for a > 0. 
This was confirmed by my simulations. Note that for 
(7 = the blockades persist and in case of random initial 
state a p oc t -1 / 4 decay can be observed |l5fl . 

The 2-PARW (a) and the 2-PARWS (b) models ex- 
hibit active steady states for a > with a continuous 
phase transition at a = 0. Therefore the exponents at 
the critical point will be those of the ARW-2 model. The 
convergence to the steady state is very slow. For a = 0.1 
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it was longer than 10 9 Monte Carlo steps (MCS). This 
has limited the simulations by approaching the critical 
point at a — 0. However as Figure |l| shows a rather good 
scaling behavior of the density versus a can be observed. 




FIG. 1. Steady state densities as the function of <r in the 
2-PARW (squares) and in the 2-PARWS (circles) models. 

The local slope analysis shows that the scaling be- 
havior extrapolates to f3 = 2.1(2) in case of the 2- 
PARWS model and to = 1.9(2) in case of the 2-PARW 
model. These values are in agreement with those of the 
2-BARW2a class (j3 = 2), where production is such that 
pair annihilation is enhanced. 
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FIG. 2. Space-time evolution from random initial state of 
the 2-PARWAS model at the critical point. Black dots corre- 
spond to A particles, red dots to B-s. 

In case of the 2-PARWAS (d) model the AB blockades 
proliferate by production events. As the consequence of 



this an active steady state appears for a > 0.3253(1) 
with a continuous phase transition. The space-time evo- 
lution from random initial state shows (Fig.||) that com- 
pact domains of alternating ..ABAB.. sequences sepa- 
rated by lonely wandering particles are formed. This is 
very similar to what was seen in case of one-component 
binary spreading processes fl33| : compact domains within 
a cloud of lonely random walkers, except that now do- 
mains are built up from alternating sequences only. This 
means that ..AAAA... and ...BBBB ... domains decay at 
this annihilation rate and particle blocking is responsi- 
ble for the compact clusters. In the language of coupled 
DP + ARW model |53| the pairs following DP process are 
the AB pairs now, which cannot decay spontaneously but 
through an annihilation process: AB + BA — ► 0. They 
interact with two types of particles executing annihilat- 
ing random walk with exclusion. 

Simulations from random initial state were run up to 
10 6 MCS. The local slopes of the particle density decay 



a eff{t) 



-\n[p{t)/ P {t/m)} 
ln(m) 



(11) 



(where m = 8 is used) at the critical point go to exponent 
a by a straight line asymptotically, while in sub (super- 
critical cases they veer down(up) respectively. 




FIG. 3. Local slopes of the density decay in the PARWAS 
model. Different curves correspond to a = 0.325, 0.3252, 
0.3253 0.3254, 0.3255, 0.326 (from bottom to top). 

At the critical point (a c — 0.3253(1)) one can estimate 
that the effective exponent tends to a = 0.19(1), which is 
higher than the the exponent of the 1+1 dimensional di- 
rected percolation 0.1595(1) Q and in fairly good agree- 
ment with that of the PCPD model in the high diffusion 
rate region (0.20(1)) @. 

In the supercritical region the steady states have been 
determined for different e = a — a c values. Following 
level-off the densities were averaged over 10 4 MCS and 
1000 samples. By looking at the effective exponent de- 
fined as 



Ine, - lne^i 



(12) 
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one can read-off: /3 e tf — » (3 ~ 0.37(2), which is again 
higher than that of the 1+1 dimensional DP value 
0.27649(4) 0, and agrees with that of the PCPD model 
in the high diffusion rate region (0.39(2)) Q . 
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FIG. 4. Effective order parameter exponent results. Linear 
extrapolation results in j3 = 0.37(2). 

Finally the survival probability (P(t)) of systems 
started from random initial condition was measured for 
sizes: L = 50,100,200,400,800,1600. The characteristic 
time t(L) to decay to P(t) = 0.9 was determined and 
shown on Figure |5[ At criticality one expects the finite 
size scaling 



t(L) oc L 2 



(13) 



where Z is the dynamical exponent. The power-law fit- 
ting resulted in Z = 1.81(2), which is far away from the 
DP value Z = 1.580740(34) (36) but close to various es- 
timates for the PCPD value Z = 1.75(10) [E>6|]2T 
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FIG. 5. Lifetime versus system size at the critical point 

In conclusion I have shown in this work that the hy- 
pothesis that I made for N-BARW models with exclu- 
sion pl[ may be extended for coupled binary production 



annihilation models. The critical point in the 2-PARW 
and 2-PARWS models, where AA and BB pairs create 
offsprings, continuous phase transition occurs at a = 
production rate therefore the on-critical exponents coin- 
cide with those of the the 2-ARW model. The simula- 
tions for the off-critical behavior of the order parameter 
have shown that the transition belongs to the 2-BARW2a 
class. The robustness of this class is striking especially 
in case of the 2-PARWS model where in principle two 
copies of PCPD models are superimposed and coupled 
with the exclusion interaction only. 

If the production is generated by different types of par- 
ticles (AB) such that alternating sequences are generated 
(2-PARWAS model) the space-time evolution will resem- 
ble to the of the PCPD model with alternating frozen 
sequences inside the compact domains. This system ex- 
hibits a continuous phase transition at a = 0.3253(1) 
with exponents in fairly good agreement with those of 
the PCPD model in the high diffusion region. In the 
model where AB pairs create offsprings in a such a way 
that prompt annihilation is possible active steady states 
are not formed for any a and the density decays with- 
out blockades for a > as p cx <~°- 5 but a crossover to 
2-ARW model scaling p oc i~ a25 occurs at a = 0. 
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